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As  a  preliminary  we  examine  a  two -ray  condition  as  shown  in  Figure  4 . 
One  wavefront  is  assumed  to  arrive  along  the  x-axis,  another  at  an  angle  0. 
The  wavefronts  are  sinusoidal  in  time  and  the  sum  at  points  along  x  is 

s(t,x)  =  A^cos(ut  -  kx)  +  A^cos^mt  -  $  -  kx  cose)  (1) 


FIGURE  4.  RAY  ARRIVALS 


The  case  treated  here  is  analogous  to  one  encountered  in  Fit  systems  with 
sinu-oidal  interference  to  a  desired  carrier. 
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A^  and  A2  represent  the  magnitude  of  the  two  wavefronts,  is  the  phase  dif¬ 
ference  between  them  at  x  •=  0,  and  k  =  2x/X  where  X  is  the  wavelength.  A 
sensor  at  point  x  responding  equally  and  linearly  to  both  wavefronts  would 
see  s(t,x).  Writing  (1)  in  envelope-angle  form, 


where 


s(t,x)  =  B(x)cos(mt  +  f  (x) ) 


B(x) 


+  A„ 


2A^2  cos  [ <f>  -  kx(l  -  cos0)] 


ij>(x)  =  tan 


-1 


A^sin  kx  +  A2sin($  +  kxcos0) 
A^cos  kx  +  A2Cos(ifi  +  kxcos8) 


(2) 

(3) 

(4) 


The  phase  obtained  using  (4)  will  be  modulo-2m.  It  is  useful  to  deal  with 
the  phase  derivative  dij>/dx  if  the  modulo-27r  ambiguity  is  to  be  avoided.  It 
can  be  shown  that 


£  =  k  -  |ci  -  cose) 


where 


1  + 


1  -  a 


1  +  a  +  2a  cos[<{>  -  kx(l  -  cos0)  ] 


a  -  A1/A2 


(5) 


(6) 


When  a.  is  large,  meaning  that  the  important  part  of  the  received  v?ave  is 
along  the  x-axis 


^  =  k  -  — „.cos-Pl  cosf,},  _  Joed  _  cos0)],  (7) 

It  fluctuates  sinusoidally  around  k  with  the  fluctuation  amplitude  decreasing 
to  zero  as  a_  goes  to  infinity.  For  a_  small 

=  k  cos8  -  ka(l  -  cosO)  cos[(J>  -  kx(l  -  cos0)],  (8) 

again  a  sinusoidal  fluctuation  which  decreases  to  zero  as  approaches  zero. 
Fot  intermediate  values  of  ji  the  fluctuation  of  dj'/dx  is  as  shown  in  Figure  5. 
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FIGURE  5 .  PHASE  DERIVATIVE  VARIATION  WITH 
DISTANCE  IN  THE  TWO-RAY  CASE. 

The  variation  is  periodic  with  period 


1  -  cos8 

and  has  peak  excursions  above  and  below  k  given  by 


(9) 


S  *  k(l  -  cos9)/ (a  -  1)  (10) 

a  =  -k(l  -  cos0)/(a  +  1)  (11) 


Figure  5  is  drawn  assuming  a  >  1.  The  fluctuation  is  around  k,  which  turns 
out  to  be  the  average  of  di^/dx.  Note  that  at  a  =  1+  8  is  positive  and  high 
in  magnitude,  and  the  fluctuation  is  highly  impulsive.  For  a  <  1,  diji/dx 
fluctuates  around  the  value  k(l-cosQ)  rather  than  around  k  and  for  a  =  1  g  is 
negative  and  high  in  magnitude,  and  the  fluctuation  is  again  impulsive  but 
negative  going.  When  is  close  to  unity  the  phase  as  a  function  of  posi¬ 
tion  (which  is  the  integral  of  dij;/dx)  is  as  shown  in  Figure  6. 

In  underwater  applications  the  angle  0  typically  found  in  long  range 
paths  is  less  than  20°.  Assuming  it  to  be  10°  the  period  is 


X 


A 

1  -  0.98 


50A 


(12) 


> 
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It  Is  interesting  to  note  that  in  analyses  found  In  the  literature  of  the  cor 
relation  distance  of  underwater  acoustic  waves,  numerical  estimates  around 
SOX  are  typically  obtained  (see  for  instance  [1]). 


t 1 ]  P.  W.  Smith,  Jr.,  "Spatial  Coherence  in  Multipath  oi  Hultimodel  Channels,” 
Journal  Acoustic  Society  of  America,  Vol.  60,  No.  2,  Aug.  1976,  pp.  305-310 


We  now  turn  to  the  more  general  case  of  n  wavefronts  arriving  at 
angles  0^,  j  «*  1,2,... .n.  The  received  sum  at  a  position  x  is 


s(t,x)  *  E  A,cos(wt  -  -  kx  cos0.) 

j=l  ■*  3  3 


(13) 


The  $  ,  j  =  l,2,...n  are  random  phase  angles  of  each  of  the  wavefronts  on 
arrival  at  the  point  x  =  0.  It  is  convenient  to  write  this  in  the  form 


s(t,x)  *  B(x)  cos[mt  +4>(x)  1 


jwt 


where 


“  Re  z(x)e 


Wvl  n  +  104  COS0J 

: (x)  -  B(x)e3vW  =  E  A  e  3  3 

j=l  3 


(14) 


(15) 


<Kx)  is  the  phase  angle  we  will  study  ar.d  as  was  done  before  we  find  the 
phase  derivative 


dip 


dx  vz(x)  dxy 


From  (15)  we  have 


dz 


~  =  ]  E  A.  k  cosO,  e  3 
dx  j-i  ^  ^ 


j(<J>..  +  kx  cosGj) 


so  that  (16)  becomes 


difr 

dx 


EA,  k  cosO.  e 


J  +  kx  cosOj) 


Im  j 


1 


J 


EA,  e 

j  3 


j  +  kx  cosOj  ) 


E  E  A  A,  k  cosO.  e 
1  j  1  3 _ 1 


jt«f.  -  ^  +  kx  (cosOj  -  cor.O^)] 


j  .  +  kx  cosO  . )  |  ? 

E  A  e  3  3 

J 


T.  E  A^A^  k  cosO^  co:;[('..  --  +  kx(cor.0j.  -  cosOj)] 

E  E  A. A  cos[iJ>  -  $  .  +  kx(cnsO,  -  ro;0.)l 

*  J  ■  J  J  J 


(16) 


(17) 


(18) 


OPK  N<>.  ?4 


MV 


6 

As  a  rule  0^  «  —  and  cose^  =  1 - j-*  Where  this  approximation  is  permis¬ 

sible  we  can  write 


where 


.  £  £  0,  a..(x) 

dx  2  11  a  (x) 

a1:J(x)  =  AjAj  cost^  -  ijij  +  kxtcosOj^  -  cos0j)] 
To  retrieve  we  must  integrate  dij>/dx;  i.e.  we  form 


(19) 


<ty(x  ) 

— j— -  dx^  =  iii(x)  -  ii<(0)  (20) 

o  X1  x 


The  integration  will  generate  the  phase  difference  between  the  phase  at  x 
and  the  phase  at  the  origin  of  integration.  From  (19)  we  see  that  one  term 
on  integration  will  te  kx,  the  linear  phase  variation  associated  with  the 
normal  phase  vs.  position  function  of  a  plane  wave  along  the  direction  of 
travel  of  the  wave.  In  beam  forming  with  an  array  of  sensors  along  x  one 
will  subtract  the  phase  progression  kx  if  the  axis  of  the  beam  is  to  be 
colinear  with  the  x  axis.  In  this  case  the  remaining  phase  difference  be¬ 
tween  a  point  x  and  the  origin  is 


-k 

2 


T.  I  0 


U 

7.  Z  a 
o  i  j 


aij(xl) 


ij(xl) 


dx. 


(21) 


2 

If  the  approximation  cosOj,  =  1  -  0^  /2  is  not  used,  the  remaining  phase 
after  correcting  for  kx  is  given  by  subtracting  k  from  (18)  and  integrat¬ 
ing  over  x. 

Numerical  evaluations  of  the  remaining  phase  difference  have  been  made 
for  a  number  of  cases. *  One  particular  case  is  shown  in  Figure  7  determined 
assuming  21  equal  amplitude  rays  arriving  at  2°  intervals  from  0  =  -20°  to 


* 

Programming  of  this  computation  and  the  one  described  later  giving  array 
pattern,  war  clone  by  Do  Juan  Ho. 
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FIGURE  7.  PHASE  DIFFERENCE  vs  POSITION 


+ 20 “  with  respect  to  the  horizontal  (see  Figure  4),  each  with  a  phase  angle 
4  randomly  selected  in  the  interval  (0,2tt).  Because  the  average  wavelength 
of  the  various  rays  as  seen  along  x  is  less  than  X,  the  wavelength  along  the 
direction  of  travel  of  the  ray,  there  is  a  linearly-tending  phase  accumula¬ 
tion  with  distance  as  seen  in  Figure  7.  On  top  of  this  accumulation  there 
is  a  random  variation.  The  fluctuation  around  a  straight  line  approximation 
to  the  phase  difference  ranges  around  +3  radians.  Thus  even  if  the  phase 
were  corrected  to  account  for  the  slope  of  the  straight  line  approximation, 
a  +  3  radian  random  error  would  still  be  encountered. 

Figure  7  was  obtained  with  one  randomly  selected  set  of  ray  arrival 
phase  angles.  Additional  examples  will  be  ultimately  computed  for  different 
sets  of  arrival  phase  angles  to  provide  data  suitable  for  obtaining  statisti¬ 
cal  averages.  Other  cases,  including  different  intervals  of  arrival  angle, 
different  ray  amplitudes,  and  different  numbers  of  arriving  arrays  will  also 
be  treated. 

Having  a  sample  function  of  phase  vs.  position, a  logical  next  step  is 
to  determine  the  gain  and  pattern  of  the  random  planar  floating  array  when 
it  is  focused  in  some  azimuthal  direction  using  conventional  beamforming, and 
when  the  source  signal  is  propagating  toward  the  array  through  the  multipath 
medium.  As  a  first  step  a  program  was  developed  for  selecting  element  posi¬ 
tions  over  a  circular  area  assuming  a  uniform  distribution  of  element  posi¬ 
tions. 

If  the  array  is  assumed  confined  to  a  circle  of  radius  P  with  uniform 
distribution  over  the  circle,  the  density  function  in  the  joint  random  vari¬ 
ables  X>  Y*  is 

t  v  1  2  2  „  2 

px  yu,y)  =  — 2  *  x  +  y  i p 

’  *  p 

=0  ,  elsewhere. 

Transforming  to  polar  coordinates,  (R,$),  we  have 


PR,*(r’*>  ~  ~~2 

ifp 
=■  0 


0  <_  r  <_  p 
0  <  <t>  <_  2  tr 

elsewhere 


QPR  No.  24 


The  marginal  densities  in  R  and  4>  are 


PR(r)  *=  ~2  0  <_  r  <_  p 

*=  0  ,  elsewhere 

P  "  ^7  *  0  <  ♦  1  2* 

*=  0  ,  elsewhere 


The  random  variables  R  and  4>  are  independent  and  independent  choices  of  these 
variables  are  made.  Sample  values  of  $  are  obtained  by  a  conventional  com¬ 
puter  program  which  selects  sample  values  uniformly  distributed  in  (0,1)  and 
multiplies  these  by  2ir.  Sample  values  of  R  are  obtained  by  picking  a  number 
Z  uniformly  distributed  in  (0,1)  and  forming 


for  then 

pR(r)  =  pz(z)  fr  =  »  0  <  r  <  p  . 

P 


Finally,  the  pairs  (r,£)  so  obtained  are  converted  back  into  rectangular  co¬ 
ordinates  by 


x  =  r  cos$ 
y  -  r  sin$ 


Using  element  positions  so  determined  the  array  pattern  was  next  found. 
The  geometry  of  the  problem  is  shown  in  Figure  8. 

Assuming  N  elements  distributed  over  the  circle,  cophased  to  form  a  beam 
along  the  y  axis,  the  array  pattern  is  given  by 

,  N  j[k(x  coutf>  +  y  sin$  -  y  )  -  a(d  )) 

*(♦>  -  |  I  I  B(dn)c  "  "  “  I 

n=l 


where  x  ,  y  ,  $,  and  d  are  defined  in  Figure  8,  and  a(d  )  is  a  phase  vs  po- 
n  ti  n  n 

sitlon  function  of  the  form  obtained  earlier  and  shown  in  Figure  7-  B(d  )  is 

n 

the  amplitude  of  the  acoustic  field  at  the.  nth  element.  ’Jhis  quantity  can  be 
obtained  using  the  earlier  analysis  but  for  our  pmposts  now  we  will  assume 


qrii  IIo.  24 


OW) 


it  constant  and  set  it  equal  to  unity  for  all  n.  Amplitude  fluctuations  as 
a  rule,  cause  minor  effects  compared  to  phase  fluctuations.  The  phase,  sample 
function  of  Figure  7,  called  now  a(x)  is  used  alone  below  to  assess  the  ef¬ 
fect  of  the  multipath  medium.  The  variable  x  in  Figure  7  is  replaced  by  d  , 

n 

with 

d  =  p  -  r  cos((b  -  $) 
n  n  n 

,  2  ,  2.1/2  ,  -1  yn  .. 

=  p  -  (x  +y)  cos  (tan - v) 

n  n  x 

n 

Computer  calculations  of  A(40,  as  described  above,  were  carried  out  for  two 

cases:  (1)  a(d  )  **  0  and  (2)  a(d  )  as  given  by  Figure  7,  and  the  results 

n  n 

are  shown  in  Figures  9  and  10.  Cate.  1  is  that  cf  propagation  through  a  trans¬ 
parent  (non-multipath)  medium  while  case  2  is  for  the  particular  multipath 
case  resulting  in  the  phase  function  discussed  above.  Note  that:  the  gains  along 
the  main  beam  in  the  two  cases  are  in  the  ratio  of  about  4.4  dB  -  a  substan¬ 
tial  factor;  the  sidelobe  structure  is  different  in  detail  but  not  in  general 
characteristics.  These  results,  it  must  be  recognized,  are  based  on  one  set 
of  random  arrival  phases  and  on  on^>  net  of  random  element  positions;  whether 
they  are  representative  remains  t<>  be  determl n<-d .  Averaging  over  many  seta 
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FIGURE  9 . 


ARRAY  PATTERN  OF  RANDOM  PLANAR  ARRAY  IN  TRANSPARENT 
MEDIUM.  ENDFIR £  BELAJ-I  FORMED  AT  AZIMUTH  ANGLE,  4  =  90°. 


RELATIVE 

VOLTAGE 

GAIN 


FIGURE  10.  ARRAY  PATTERN  OP  RANDOM  PLANAR  ARRAY  ASSUMING  MULTIPATH 
PROPAGATION.  ENDFIRR  REAM  FORMED  AT  AZIMUTH  ANGLE  $  *  I 
FOR  TRANSPARENT  MEDIUM. 

of  arrival  phases  and  element  positions,  as  well  as  carrying  out  additional 
computations  with  other  system  parameters,  remain  to  be  done. 

Fred  Haber 
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